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Abstract: Standard perturbation theory (SPT) for large-scale matter inhomogeneities is 
unsatisfactory for at least three reasons: there is no clear expansion parameter since the 
density contrast is not small on all scales; it does not fully account for deviations at large 
scales from a perfect pressureless fluid induced by short-scale non-linearities; for generic initial 
conditions, loop corrections are UV-divergent, making predictions cutoff dependent and hence 
unphysical. 

The Effective Field Theory of Large Scale Structures successfully addresses all three issues. 
Here we focus on the third one and show explicitly that the terms induced by integrating 
out short scales, neglected in SPT, have exactly the right scale dependence to cancel all UV- 
divergences at one loop, and this should hold at all loops. A particularly clear example is 
an Einstein deSitter universe with no-scale initial conditions Pi n ~ k n . After renormalizing 
the theory, we use self-similarity to derive a very simple result for the final power spectrum 
for any n, excluding two-loop corrections and higher. We show how the relative importance 
of different corrections depend on n. For n ~ —1.5, relevant for our universe, pressure and 
dissipative corrections are more important than the two-loop corrections. 
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1. Introduction 

On very large scales, inhomogeneities in our universe are small, of the order 10 -5 . On small 
scales, due to gravitational collapse, clumps of matter form. Numerical methods such as 
N-body simulations can tackle the full non-linear problem and therefore have been widely 
employed to describe all scales. In recent years, with the development of the Baryon Acoustic 
Oscillation (BAO) technique as a probe of the expansion history of the Universe, interest has 
grown in the so-called mildly non-linear regime of structure formation. These are scales where 
non-linearities are small and thus amenable to perturbative methods but still non-negligible 
to interpret precision measurements. Numerical simulations can be used to study this regime, 
although to probe cosmological models whose parameter space has a large dimensionality re- 
quires significant computational resources. In addition, numerical simulations may tend to 
make the underlying physics opaque, or at least one is much more satisfied when one can 
achieve a relatively accurate analytic understanding of what the simulations show. This is 
why, despite the impressive development in the computational speed, there is still a large ef- 
forts to improve analytical methods and thus also improve our understanding of the evolution 
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of structures. Most analytical methods are based on a perturbative expansion in small inho- 
mogeneities and therefore are valid only in between the Hubble scale H and some non-linear 
scale fc/vx ^ clH, with a the FLRW scale factor. Much work, including the present one, has 
been devoted to pushing the accuracy of the analytical methods closer to the non-linear scale 
and understanding its physical origin. 

In the standard (Euclidean) perturbation theory (SPT) (see [1] are references therein) the 
final results, e.g. for the density power spectrum, are organized in a loop expansion where the 
small parameter is the initial power spectrum of inhomogeneities. There are at least two con- 
ceptual issues worth emphasizing: first inhomogeneities are large at small scales, which makes 
it hard to rigorously define the expansion parameter. Furthermore at sufficiently small scales 
multiple streams develop so the fluid equations typically used become invalid. Second, for 
generic initial conditions, loop integrals are typically divergent, but no counterterms appear 
in the theory with the appropriate scaling to cancel them. This means that the predictions 
depend on the cutoff and are hence unphysical. The sensitivity to the cutoff also points to 
the fact that the answer depends on small scales where the physics in not being modeled 
properly. 

Both issues are successfully addressed within the Effective Field Theory (EFT) of Large 
Scale Structures (LLS) [2, 3, 4], henceforth EFToLSS. The main idea is the usual underlaying 
principle of EFT. Integrate out the short-scale modes and systematically incorporate their 
effect on the large-scale dynamics (see also [5] for a similar approach). Starting with a 
configuration of collisionless dark matter particles in the Newtonian limit in an expanding 
universe, one derives a system of fluid equations. The effects of the short scales is to induce 
a non-vanishing speed of sound, dissipative corrections and a stochastic noise on top of the 
perfect fluid that one would have obtained if these effects had been neglected. On the one 
side, this is remarkable if one remembers that we started with collisionless particles. On the 
other this is precisely what we should have expected based on the EFT philosophy: all terms 
compatible with the symmetries of the problem should be present. These additional terms 
encode the effects from the small scales whose dynamics is not captured by the EFT on the 
large scales which are now correctly modeled. 

The goal of this paper is to clarify how the approach of the EFToLSS solves the second 
problem mentioned above, namely the fact that SPT predictions are cutoff dependent and 
hence unphysical. Building on the results of [3, 4], we show in detail that the effective 
terms induced by integrating out the short modes have exactly the right scale and time 
dependence to cancel the divergent terms. This result is general and makes a strong case 
that the EFToLSS, rather than SPT, is the theoretically consistent way to perform Euclidean 
perturbation theory. In order to show this in the clearest possible way, starting in section 2, 
we focus on an Einstein de Sitter (EdS) cosmology containing only non-relativistic matter. 
This has a few advantages: first EdS is a good description of our universe after matter 
radiation equality with the largest deviations arising at late times when Dark Energy starts 
dominating; second, thanks to the simple form of the growth factor in perturbation theory 
almost all results are amenable to analytical computations; finally, provided the right initial 
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conditions are specified, namely a no-scale initial power spectrum P$ ~ k n , the system enjoys 
a self-similar evolution, which ensures that analytical solutions can be constructed that are 
valid at all times. 

We discuss the renormalization of the theory in section 3. We employ both a hard cutoff, 
which breaks self-similarity, and dimensional regularization which preserves it. The final 
results agree with each other and satisfy self-similar scaling as expected. For EdS we write 
down explicit formulae for the final non-linear density power spectrum (without including 
two-loop terms and higher) including the leading effective corrections for any n. Our final 
result is given in (3.6). We also highlight the relative importance of the different terms as the 
non-linear scale is approached (see figure 1). Our expression (3.6) contains zero, one or two 
free parameters (excluding two- loop corrections or higher), depending on n. In the specific 
case n ~ —1.5, which is appropriate for our universe close to the present non-linear scale, only 
one parameter appears up to two-loop correction, namely the one associated with the speed of 
sound (whose phenomenological effects were studied in [3, 4]). We present some preliminary 
comparison with results from published numerical simulations in EdS for n = —1.5 and 
n = — 1 in section 4 and show that there is good agreement, but clearly more work is required 
to make a more convincing comparison. We conclude in section 5 with a discussion and ideas 
for future research. 



2. Smoothed perturbation theory 

Accounting for the effects of integrating out small scales, the equations of motion for the 
density contrast 5 = p/p — 1 with p the average density, the velocity v % and the gravitational 
potential </> in the Newtonian limit (neglecting vorticity, which decays at linear order) around 
an FLRW universe are [2, 3] 

a r s l + d i [(i + Si)vj\=o, (2.i) 

Q r 2 4 2 i 2 

d T v\ + Hvj + di4> + vjfdkv} = -c 2 s d i d l + ~~d 2 vi + bv Z, av d i d j v 3 l - A J 1 ... , (2.2) 

4 rl Art 

did^l = = = 3 J3^ 6l (2 .3) 

2 2 2 a 

where r is conformal time, H = d T a/a = aH with H = dta/a and Q. m (r) = pj'ili 2 M 2 - v with 
Slo = f2(ao = 1). Some additional comments are in order. The subscript I indicates that 
the corresponding field has no short-wavelength fluctuations, i.e. it has been smoothed on a 
certain scale A. These equations can be solved perturbatively when A < h^L (we will define 
kjvL in (2.25)), so that the smoothed long-scale perturbations are never large. On the right 
hand side of (2.2) one finds the corrections induced by the short scales. This is a double 
expansion in perturbations and derivatives acting on them. At zeroth and linear order in 
perturbations to lowest order in derivatives we have four terms. Each one of them comes 
with a infinite series of higher derivative corrections which we include in the ellipsis. c 2 ,c 2 v 
and c 2 v are real time-dependent functions that physically correspond to the speed of sound, 
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bulk and sheer viscosity respectively. AJ l is a time and scale dependent random variable 
representing stochastic noise induced by the fact that any specific short-scale configuration 
deviates somewhat from the ensemble average. The values (and distribution function in the 
case of AJ l ) of these parameters are not predicted within the EFToLSS approach. Rather, 
they can be computed in the microscopic theory or, when this is too hard, they should be 
extracted by comparison with observations or numerical simulations. 

2.1 Loop corrections in an Einstein de Sitter universe 

One of the main goals of the present paper is to show that the new terms induced by integrat- 
ing out the short scales have exactly the right scale dependence to cancel the UV divergences 
arising in perturbation theory. Because the EFToLSS is derived by including all possible 
terms compatible with the symmetries of the problem (homogeneity and isotropy of the back- 
ground, conservation of mass, etc.) this should be true in any cosmology at all loops. For 
concreteness we consider an EdS universe because the formulae simplify considerably and 
because of its relevance for our universe. We also neglect two- loop and higher terms, but 
expect that our results generalize straightforwardly to higher orders. 

Consider a flat universe filled only with non-relativistic matter, a.k.a. an Einstein-deSitter 
(EdS) universe. Then VL m = 1 and one finds 

% = -, a=(-) , n = nVa~, n' = ~. (2.4) 
t \t J 2a 

Neglecting for the moment the effective terms in the right-hand side of (2.2), the system 
reduces to the usual SPT equations whose perturbative solutions are well known (see e.g. [1]). 
Here we simply quote the one-loop result for the correction to the power spectrum 

P(a, k) = P lin (k, a) + P 22 (k, a) + P 13 (k, a) + P^(k, a) + Pj(k, 



with 



,a), (2.5) 

Ptin = a 2 P in , (2.6) 



a 4 k 3 



Pis(k) = ^^2 p m(k) I drP in (kr) 

(12 3 
- 158 + 100r 2 - 42r 4 + -3 (r 2 - l) 3 (7r 2 + 2) In 



1 + r 



1 — r 



(2.7) 



P22ik) = 984^ J dr J dX Pm{k r ) P ^ 1 + r ~ 2ra ) (1 + r 2 _ 2rx) 2 ■ (2-8) 

Before discussing the new lowest order corrections induced by the short modes, namely 
P c 2 and Pj, let us pause and comment on the one- loop terms Pi3,22- These can be both 
UV and IR divergent depending on the initial power spectrum under consideration. For no- 
scale initial conditions Pj n = Ak n with A a constant, the conditions for the occurrence of 
divergences are summarized in table 1. For the moment we will ignore these divergencies but 
will discuss them in detail in section 3. 
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To compute these corrections we notice that at linear order in perturbations one can find 
a second order differential equation for Si (valid in any FLRW cosmology) 

- a^dfo - a (2H 2 + aHH >a ) 8 a 5 l + ~ft 2 ft m ^ = -c 2 s d 2 5i - c 2 v ad 2 d a 5 l - J , (2.9) 

where we used the scale factor a as time variable and defined c 2 = c 2 v + c 2 v and J = d{AJ l . 
Let us treat the terms on the right-hand side perturbatively. Since the left-hand side of 
equation has no spatial derivatives, the zeroth-order solutions is given in the form 

S t = 5i(k) [C+D + (a) + C-D-(a)] , (2.10) 

with C± two integration constants and 5i(k) a stochastic variable fixed by some boundary 
condition. The growing and decaying solutions are D + = a and D- = a -3 / 2 and in the 
following we will discard the latter. The retarded Green's function in EdS (we discuss the 
d-dimensional generalization of this result in appendix B) takes the form 



G(a,a) = 9(a-d)ln 2 

Notice that 



, 3/2 



I 



a a 



a n+l 



(2.11) 



daG^ = - nln{n + b/2y (2.12) 
which explains the well know structure of SPT in EdS 5 n oc a n . 
2.2 Self-similar solutions 

In an Einstein de-Sitter cosmology the equations of motion posses a scaling symmetry. For 
simplicity let us consider the case when matter can be treated as a cold fluid. Although 
this is not the case on small scales and we should be studying the dynamics of collisionless 
particles, the equations for a collection of non-interacting particles has the same symmetry 
we now discuss. The equations are just those in (2.2) neglecting for the moment the effective 
terms in the right hand side. 

d T 5 + di [(1 + 6)1?] =0, (2.13) 
d T v i + Hv i + d i (l> + v k d k v i = 0, (2.14) 

W = ^H 2 S, (2.15) 

where H = 2/r. One can take a solution of these equations [6(x, t),v 1 (x, t), cp(x, r)] and scale 
it to get another one: 

S(x,t) = 5(\ x x, X t t) , (2.16) 
v i (x,T) = ^v i (X x x,X r T), (2.17) 

Ax 

0(x,r) = (y) 4>{X x x,X T r) . (2.18) 
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That this is a solution can be checked by simply replacing it in eq. (2.13). Although in this 
way we obtain a new solution for any choice of A^ and A r , not all the choices correspond 
to a realization of the same cosmology. By this we mean that not all solutions have initial 
conditions sampled from the same power spectrum. By making a specific choice of the relative 
x and r scalings one can get a family of solutions such that for any A they are just different 
realizations with the same initial statistical properties. 
Consider A 2 (/c,r), 

A 2 {k,r) = fc3p ( fc ' T ) t with (5(k)5{k')} = {27rf5 3 (k-k')P(k,T). (2.19) 
After scaling we have: 

A 2 (k, t) = A 2 (k/X x , A r r). (2.20) 

We will first study this relation at very early times such that all quantities can be com- 
puted in linear theory. We will assume that [5(x,t),v*(x,t),^>(x,t)] were sampled from a 
Gaussian distribution. The time dependence of the scaled solutions is the same as the original 
one so both are in the correct growing mode solution. The only remaining issue is to make 
sure that the scaled solutions have the same power spectrum at early times. Assuming power 
law initial conditions for 5 with index n, i.e. A 2 niUal (k,T) oc r 4 fc n+3 , we get 

A 4 

^■initiai{k-, T ) oc n + 3 A initial (k, t). (2.21) 

Thus if we pick 



4 

Ax = A" , (2.22) 

then the different rescaled solutions are just different samples of the same initial power spec- 
trum. If this is so, then at any given time the power spectra of all the scaled solutions (which 
is an ensamble averaged quantity) should be the same as the original un-scaled solution. 
Thus, 

A 2 {k,r) = A 2 (k,r) = A 2 (k/X^ , Ar). (2.23) 

This is satisfied if 

A 2 (k,r) = A 2 (k/k NL ), (2.24) 



where 



2tt 2 



k NL = 1^2 « ^" 4 > ( 2 - 25 ) 



Aa 

so that k/k^L is invariant under the scaling (2.22). With this definition, Ajj n (fcjvx) = 1- 
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Perhaps this is just a complicated way of saying that the only scale in the problem is 
the non- linear scale and thus the answer should just be expressed in terms of it. Of course 
in practical situations, such as for example when running a numerical simulation with finite 
computer resources, one is forced to introduce additional scales such as the size of the box 
or the inter-particle separation. As long as these scales are sufficiently removed from the 
non- linear scale their effects should be small. This has been checked in practice for a wide 
range of re's. Convergence is progressively more difficult to achieve the closer the spectrum is 
to n = —3 at which point all scales become non-linear at the same time. 

Given that the full solution satisfies the scaling symmetry, both at long and short scales, 
the new terms appearing in the smoothed equations of the EFToLSS (2.2) and (2.9) should 
also satisfy it, since they arise from integrating out the short scales. This implies that each 
term should scale according to its dimensions. For example: 

c 2 ,c 2 bv ,c 2 v oc (2.26) 

T 

Since these terms represent averages over the short wavelengths, the scaling symmetry should 
only be respected if we consider the combined x — r rescaling that leaves the initial spectrum 
fixed, namely (2.22). We can then use (2.22) and (2.9) to deduce the appropriate time 
dependence of the sound speed coefficients, 

c?,cL,cLocrW. (2.27) 

This scaling of c 2 can be derived as well by requiring that the Jeans scale is a constant fraction 
of the non-linear scale, since again there should be only one independent scale in the problem. 

The properties of J in (2.9) can be similarly obtained. In this case the power spectrum 
of J should scale in such a way as to preserve the scaling symmetry under the combined x — r 
scaling. The only additional complication now is that J, unlike c 2 , is spacetime dependent, 
and hence can contain and unknown function of h/h^L, which, being invariant under the 
combined x — r scaling, can not be fixed with the arguments we have used so far. In particular, 
by inspection of (2.2), one finds 

A 2 j(k,T) = T-*f(k/k NL ), (2.28) 

where / is some arbitrary function. In the limit of A; <C fc/yx the k dependence can be 
constrained following arguments in [10]. Mass and momentum conservation imply that on 
large scales the corrections to the density power spectrum from these terms scale as k A rather 
than k° or equivalently: 

A 2 j(k, r) a T~ 4 (k/k NL ) 7 (k « k NL ). (2.29) 
2.3 The self-similar power spectrum 

Now that we have the scaling of the effective coefficients, we can derive a formula for the power 
spectrum. For this section we will ignore the issue of divergences, postponing a thorough 



- 7 - 



discussion to the next section. Because of self-similarity, we expect that the final A 2 ot (and 
not Ptot) should be a function of just k/k^L- For this to be true the k dependence of the 
various terms will have to combine with the r dependence in a precise way. This will provide 
us with a cross check of the scalings computed in the previous section. 

Let us start considering the corrections coming from c 2 and c 2 . Their leading order 
contribution to 5 (from now on we drop the label I for large scale perturbations, leaving it 
implied that all the fields we will be discussing have been smoothed above a certain scale 
A < k^i) perturbations is given by 

5 C 2 = J daG(a, a)k 2 (c 2 + c 2 v ad a ) 61 , (2.30) 

where 5\ is the linear order perfect fluid solution, G the Green's function in (2.11) and in EdS 
a = (t/tq) 2 . The linear contribution of this term to the power spectrum is hence found to be 

A 6 2 2 oc k 5 (c 2 + c 2 v )P in (k)r 4+2 , (2.31) 

where the extra r 2 comes from the Green's function. For self-similar initial conditions Pi n = 
Ak 11 , the exponent of k is 5 + n, which perfectly combines with the exponent of r, once 
the scaling (2.27) has been taken into account, to give (k/k^L) 5+n , up to higher derivative 
corrections. 

For the noise term things work out analogously. At linear order one has 



daG{a,ti)J, (2.32) 
and hence 

A^ocA^AOr 4 , (2.33) 

where the r 4 ~ a 2 comes from the Green's function. This again combines perfectly to give 
(fc/fcjvx) 7 , up to higher derivative corrections. Finally, let us discuss higher loop correc- 
tions. Since loops count the exponent of the linear power spectrum, one expect A 2 v _ i cx 

(fc/W n+3)(iV+1) - 

In summary, the self-similar power spectrum is given by 

AL = AL + ^N-ioop + A^2 + A 2 loise , (2.34) 

N 

AL= (t^-T + \ (2-35) 



& 2 n-ioo P « ( irH . ( 2 - 36 ) 



kNL 
k 



NL 
k 



(n+3)(iV+l) 



n+5 



A c 2 oc ( - — J + higher derivatives (2.37) 
k \ 7 

A-noise x I i ) + higher derivatives . (2.38) 



kNL 

and it is just be a function of k/kNL as anticipated. 
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3. Cancellation of divergences and final power spectrum 



It is time to confront the issue we neglected so far. As we stated in table 1, the total one- loop 
contribution is UV divergent for n > — 1 and IR divergent for n < —3. For n > —3 both P\% 
and P22 are IR divergent but their divergences cancel each other exactly. This is expected to 
happen at all orders and was discussed more in details in [11]. In this paper we simply focus on 
n > —3 and safely ignore the issue of IR divergences 1 . UV divergences instead appear naively 
more problematic (see e.g. [6]). First, in order to get finite results one needs to regularize the 
theory. The details of how things work out depend on the regularization procedure and we 
will discuss cut-off regularization and dimensional regularization separately. 
The final result without including two- loop terms and higher is: 



3+n [ / 7 \ 3+n 



k y +n ( k 



A z = ^ <l + 



a(n) + a(n) In 



k 




kNL ) \ V k NL 

( k \ 5+n 
+/^( n ) I ~, ) [1 + higher deriv.] + jOn) 

\kNLj 

where /3 and 7 are fitting coefficients, a and a can be computed in SPT (see table 3.1) 
and higher derivatives come with even positive powers of k/k^i. We decompose 2 a(n) = 
&22(n) + a\z{n) and the same for 5, with 

0*13,22 + «13, 22 In ^— j = » ( 3 - 2 ) 

and give the explicit values in table 3.1. Note that, because we did not compute two- loop 
and higher terms, only terms that scale with a lower power of k/k^L than the two-loop term, 
namely (k/kNL) 3 ^ n+3 \ should be included in a consistent calculation. Figure 1 shows the 
k/k^L scaling of the various terms as a function of n. For example in the absence of two loop 
terms it is only consistent to include noise corrections for n > —0.5 and the sound speed terms 
for n > —2. For n relevant to the non-linear scale in our late universe, i.e. n ~ —1.5, the sound 
speed corrections are more important than the two-loop terms but the noise corrections are 
even less important than three loop terms. Finally notice that for those n's for which log's 
appear, the finite loop term has the same ^-dependence as either the noise or the speed of 
sound corrections (we will discuss this more in the following). In these cases a is degenerate 
with either 7 or /3, which are fitting parameters, and hence its value is irrelevant and we omit 
it. 



1 For the case n = — 1 there are IR divergences that leak into the computation of Pi3,22 when performed 
using dimensional regularization. As all other IR divergences for n > —3 also these log-divergences cancel 
exactly once we add up the two one-loop contributions. 

2 This is related to the as in [1] by Qhere = Q-s CGS F [(3 + n)/2] |. For this conversion it is important 
to notice that eq. (168) of [1] is not valid in general, but only for n = —2. The general relation is instead 
(k NL R y i+n = r [(3 + n)/2] § 
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UV div 


IR div 


Pl3 


n > -1 


n < -1 


P22 


n > 1/2 


n < -1 


Ptotal 


n > -1 


n < -3 



Table 1: The table shows for which n UV or IR divergences arise at one- loop order. 



3.1 Cut-off regularization 

Let us now discuss how to regularize the theory by introducing a UV cutoff on the momenta. 
This choice is very natural in the EFToLSS approach since the fields have been smoothed over 
scales of order A from the beginning. Then all loop integrands, which include positive powers 
of the power spectrum, go quickly to zero when the argument of any of the power spectra is 
much larger than A. This makes the loop integrals finite and cutoff dependent. This cutoff 
is time independent in the derivation of EFToLSS and hence breaks self-similarity explicitly. 
This is not a problem per se, since as we will see the final results after renormalization are 
cutoff independent and hence preserve self-similarity as expected. In fact a regulator that 
does not respect the symmetries of a theory is not uncommon in (quantum) field theory; it 
is a consistent procedure but it makes the steps preceding the final result less transparent. 
In the next subsection we consider dimensional regularization, which preserves self-similarity 
throughout. 

The one-loop divergences can be seen by expanding (2.7) and (2.8) for large q. Up to 
numerical coefficients one finds 

P 13 ~ a 4 k 2 P n (k) I dqPn(q), (3.3) 
P22 * J* ^Pn(q) 2 , (3.4) 

from which the results for the UV divergences in table 1 can be easily extracted. These two 
terms have precisely the same k dependence as P c 2 in (2.31) (as already discussed in [3, 4]) 

and Pnoise m 

(2.29) and (2.33), respectively (A 2 was defined in (2.19)). This means that any 
A dependence arising in the two one-loop integrals above can be absorbed by a counterterm 
in the form of a specific c 2 (A) or Aj(A). This is a crucial point so let us stress it again. The 
additional terms present in the EFToLSS, which arise from integrating out the short scale 
degrees of freedom and which are neglected in SPT, are precisely the terms needed to cancel 
the cutoff dependence and hence make physically meaningful predictions (this was recognized 
already in [2, 3, 4]). 

It can be easily checked that the counterterms needed to cancel the divergences do not 
have the correct time dependence required by self-similarity, but this is to be expected. Since 
the cutoff regularization breaks self-similarity explicitly by introducing another independent 
scale A, also the counterterms violate it. The important point is that, after the counterterms 
have canceled the divergences, the final results respect self-similarity. 
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Figure 1: We show the exponent of fc/fcjvL for various contributions to the power spectrum as 
function of the power in the initial conditions Pu n oc k n . The full orange lines are zero, one and two 
loop contributions from SPT (from bottom to top); the black dot-dashed line is the term A c 2 while 
the dotted red line refers to Aj. The three background regions labeled 1-loop, c 2 8 and J indicate which 
is the most important correction to An n at that value of n. For our universe, near the non-linear scale 
-1.5 (eg. [7]). 



After renormalization the only terms that survive are the finite cutoff independent parts 
of the loop-integrals, which can be computed using the expressions (2.7) and (2.8) (or (A. 9) 
and (A. 10) using dimensional regularization), and the finite parts of (? s and Aj, which scale 
as discussed in subsection 2.3. Since the details of the computation have been presented 
elsewhere [6] and are lengthy and not particularly illuminating, we do not present them here, 
but we attach to our publication a mathematica file [12] that computes, as an example, the 
finite part of P13 for —3 < n < 4 (P22 is slightly more difficult because of the angular integral, 
but can be computed with some work [6]). 

The last point to discuss concerns logarithms. For some n's, the loop integrals have a 
logarithmic divergence, or equivalently a term ln(A/fc). As customary in field theory, the 
In A is canceled by a couterterm, but the lnfe remains in the final result. These log terms at 
one loop come with the same k power law as the finite terms, i.e. {k/kNL) 2 ^ +n \ as expected 
by dimensional analysis. This can also be seen by looking at the scale dependence of the 
most divergent terms in (3.3) and (3.4), i.e. k 5+n A n+1 and /c 7 A 2n_1 (we added a factor of 
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-3/2 
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-1/2 
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3/2 


2 


5/2 
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a 13 


112 


992-7T 
6,615 




416-rr 
8,085 


2 

TV 

336 








2 

7T 

168 






«22 


75-it- 1 
784 


-0.232 




.698 


29-7T 2 

784 








7T 

392 






«13 








61 
315 











4 

105 











20 
1,323 


"22 

















9 

98 





31 
16,464 





359 
26,880 





a 


1.38 


.239 




.537 


.336 








-.0336 






a 








.194 








-.0918 


.0381 


-.00188 





-.0134 


.0151 



Table 2: In the table we report the coefficients of the finite terms in the one-loop corrections to the 
power spectrum, computed in perturbation theory with dimensional regularization. We omit those 
cases in which a is degenerate with a fitting parameter. We attach to this publication a Mathematica 
notebook [12] that generates these numbers and allows to compute their analog for any n. 



k 3 as appropriate for discussing A 2 rather than P). When the subleading terms scaling 
as positive powers of (k/A) are included, A becomes k and one obtains again the scaling 
claimed above. This means that for the counterterms to be able to cancel the In A in the log 
divergences, one needs A 2 2 or A 2 - or their higher derivative corrections to have exactly the 
same ^-dependence as the loop corrections. For the one-loop corrections considered here, this 
happens when 2(n + 3) = 5 + n + 2m or 2 (re + 3) = 7 + 2rre for some positive integer m, since 
higher derivative corrections come with an even number of additional powers of k [2] . Hence 
we expect one-loop log divergences for n = — 1 + 2m and n = 1/2 + m, with m a positive 
integer, from P13 and P22 respectively. This can be confirmed by direct computation 3 [12] 
and can be seen in table 3.1, where ayzpi are non-vanishing exactly for re = —1 + 2m and 
re = 1/2 + m, respectively 4 . 

So, after subtracting the divergences with appropriate counterterms, the loop corrections 
take the form 

„ 2 / k V+" / k y (s+ "> r , , „ / k \ 

Here a and a can be computed directly and are given in table 3.1. The direct computation 

3 Notice that using the formulae in appendix B of [6] one finds a log divergence in P22 for n = 1, which 
would be in contradiction with the discussion above. The reason is that there is a typo in [6] in the sign of the 
fourth term in P22, which should be —11/4 instead of plus. We re-derive this result in appendix A.l. After 
correcting the typo the log divergent n's from dimensional regularization perfectly match the expectations 
from the cutoff regularization discussed here. 

4 One subtlety arises for n = —1 (see also footnote 1), when there are logarithmic IR divergences, which 
contribute to Pi.3,22 computed in dimensional regularization as discussed in appendix A. For this case one 
can check that (S22 = 1/3 = —013 are exactly the expected IR logarithmic divergences, e.g. by computing the 
loop integral with an IR cutoff. E.g. one can see this in table 5 of [1], keeping in mind that our numbers are 
a factor of (2tt) 3 / (2ty 2 tv) = 4 lower than theirs. The {2n) 3 comes from the different Fourier conventions, the 
factor 2-7T 2 from the fact that we present our results with respect to fcjvz, and finally the last ir is because of 
the normalization chosen in that table P/(nA 2 a 4 ). In table 3.1 we have subtracted these log-divergences from 
the values computed in dimensional regularization 0:13 = —44/315 and 022 = 1/3 



(3.5) 
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using (2.7) and (2.8) is in principle straightforward but rather lengthy, especially for ^22- In 
the next subsection we will show how, using the results of [6], the a's can be easily computed. 
(3.5) was given already in [6] for n < —1. The result of our approach shows two qualitative 
differences. First there are additional corrections to the power spectrum, that, summing up 
(3.5), (2.31) and (2.29), takes the form 




k x3+n 



kNL 



a(n) + d(n) In 



k 



kNL 



k \ 5+n / k x 7 



(3.6) 



+/^( n ) 1 [1 + higher deriv.] + 7(71) [1 + higher deri v.] . 

\ kNL J \ kNL J 

Second this result is valid for any n > —3 (below which IR-divergences become problematic) 
and not just for — 3 < n < — 1. No divergences appear in the final result and self-similarity is 
preserved for any n > — 3. 

Before concluding, let us remark that the InkNL term in (3.5) comes from A 2 2 or A? ; ■ 
when choosing a particular time dependence (since kNL is a function of time) of c 2 . and c 2 or 
Aj. As anticipated, this time dependence is different from the one we deduced in subsection 
2.2 using self-similarity. The reason is again that since we are using a regularization that 
violates self-similarity, so do the counterterms. Notice that the final result 3.5 is self-similar 
as it should. 



3.2 Dimensional regularization 

Dimensional regularization applied to SPT was first discussed in [6]. There, formulae were 
provided for the finite and log-divergent one-loop corrections in three dimensions, since the 
power law divergences disappear in dimensional regularization. We re-derive those results 
in any dimension (correcting a typo in P22; see footnote 3) in appendix A and we have 
included them in an attached mathematica file [12] for the convenience of the reader. One 
advantage of this regularization procedure is that it respects self-similarity, so that all terms 
obey the scalings that we computed in section 2.2. Using the dimensionally regularized one- 
loop integrals ^22,13 in appendix A.l, we have computed the numerical coefficients of the 
finite one-loop corrections for the n's shown in table 3.1. Notice that for each n only one such 
coefficient is needed in the final formula, since when the log appears, i.e. when d ^ 0, then 
some fitting parameter, either /3 or 7, is degenerate with a, whose precise value we omit since 
it is irrelevant for the final shape of the power spectrum. 

Logarithmic divergences arise for those n's for which (A. 9) and (A. 10) are divergent. This 
happens for n = —l + 2m and n = 1/2 + m, with m a positive integer, in agreement with what 
we found in the previous subsection (see footnote 3). In these cases, one needs to expand 
the (f-dimensional version of these formulae around the dimension d = 3 by substituting 
d — > d + e and sending e to zero. Unfortunately [6] gives only the results for d = 3, rather 
than as function of d. We re-derived the general formulae in appendix A.l, the result is given 
in (A. 7) and (A. 8). Schematically, the leading order of a typical term in (A. 7) or (A. 8) looks 
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like 



r(...)r(...)r(...) / k 



2n+3+e 



r(...)r(...)r(...) \k NL 



(3.7) 



where the argument of the T's depend on d and n (for P13 there are only two T's but this 
is irrelevant for the following argument). The fact that this expression can be written as 
function of k/k^L for any d (and n) is a consequence of self-similarity, which, as we check 
in appendix B, holds in any number of dimensions (with an appropriate d- and n-dependent 
definition of k^L as in (B.13)). 

Logarithms arise when the argument of the T's in the numerator (since T never vanishes, 
no divergences come from the denominator) become a null or negative integer as d — > 3. Then 
the expression has a pole 1/e with e = d — 3. Using 



where # are the specific numerical factors which determine the precise n-dependent value of 
a (see table 3.1). Notice that the finite terms in (3.9) (determining the coefficient a of the 
finite loop contribution) are always degenerate with some fitting parameter (/? or 7 in (3.6)) 
and are therefore irrelevant for the shape of the power spectrum. This is a consequence of the 
fact that log-divergences arise only when the loop corrections have the same k dependence 
(hence the degeneracy) as the noise or speed of sound corrections, as explained in the previous 
subsection. Notice that, since the In kjyL terms come already from the one- loop integrals, the 
effective coefficients in the terms A^ 2 and A^ oise now have precisely the time dependence that 
we derived in subsection 2.2 using self-similarity. This was to be expected since dimensional 
regularization respect the scaling symmetry (see appendix B). 

4. Comparison with numerical results 

In this section, we compare the final predictions of the renormalized EFToLSS with numerical 
simulations, for self-similar initial conditions in EdS, (3.6). 

The first point to make is that the EFToLSS makes a definite prediction for any self- 
similar initial power spectrum Pi n {k) = Ak n with any n, as one expects from a self-consistent 
theory. On the contrary, SPT provides a prediction only for — 3 < n < — 1, which highlights 
some inconsistency in the theory. This is a key point, independent of any phenomenological 
discussion, so let us stress it again. The EFToLSS, as opposed to SPT, is the theoretically 
consistend way to study cosmological perturbation theory because it makes physically mean- 
ingful predictions (i.e. cutoff independent) for any initial condition 5 . Second, as it is common 

5 Strictly speaking here we only discuss self-similar initial conditions in EdS, but the fact that the EFToLSS 
has the right counterterms to cancel loop divergences is general. 




+ 0(e) 




with m a non-negative integer, one can expand (3.7) as 




(3.9) 
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in effective field theories, the predictions of the EFToLSS have fitting parameters, correspond- 
ing to the renormalization conditions that needs to be imposed to select a specific solution of 
the renormalization group flow. As a consequence, the EFToLSS predictions will always fit 
the numerical (or observational) power spectra as well or better than SPT. To study whether 
the improvement is just due to the additional fitting coefficients or it really comes from cap- 
turing the relevant physics requires a dedicated analysis that accounts for the error bar in 
the simulation (observations). In addition, one might not be particularly impressed by fitting 
a smooth line such as the power spectrum with a polynomial. More convincing evidence 
would come from a realization by realization comparison and from a direct measurement of 
the short mode stress tensor and thus the EFToLSS parameters in simulations. This has 
been done in [3] for the ACDM cosmology. We leave this for future work. Finally, when the 
coefficients (3 and 7 are fixed by fitting the power spectrum, there is some arbitrariness in how 
the fitting is performed. We choose to perform a least x 2 fit over the range .01 < k/k^L < 2 
approximating the error bars in the simulation to be /c-independent. We chose this particular 
range according to the following considerations. The lower bound comes from the size of the 
simulations. The upper bound should be of order one, but there is some ambiguity in how 
the non-linear scale is defined. We used that A 2 (/ctvl) = 1 but other choices are possible, 
which differ from this by some order-one factor. Operationally, one would like to define fc/vx 
by computing higher and higher order corrections and looking for the k/k^i beyond which 
the fitting cannot be improved. As a rough estimate of this point we take k/k^L = 2. 

Now that we have made these cautionary remarks, we can move on and discuss how 
(3.6) compares with simulations. We use the fitting formulae given in equations A2 of [8] 
(for n = —1.5,-1) and in 22 of [9] (for n = — 1). It would be nice in the future to consider 
other cases as well. For each n, we plot in figure 2 the fractional difference of EFToLSS (red, 
thick line) and, when available, SPT (blue, dashed line) predictions from the fitting formulae 
extracted from the simulations. In the left and right panels we show the same functions but 
with a linear or log abscissa, respectively. As shown in figure 1, for —3 < n < — 1 the leading 
correction to the linear power spectrum comes from the one loop term. All of these are known 
and summarized in table 3.1. For —2 < n < — 1 the next-to-leading comes from A 2 2 . For 
n = —2 the A 2 2 correction is of the same order as the two-loop term, and, as discussed in the 
previous section, we would expect a log correction to arise (at two loops, i.e. multiplied by 
(fc/fcjVL) 3 ^ 3+n ^ =3 )- Since the coefficient of the finite two-loop term are not known and hard to 
compute we cannot show plots for this case, even though the results from simulations would 
be available 6 [9]. Instead for n > —2 the two- loop correction is at least the next-to- next-to- 
leading order and therefore negligible. Notice that n ~ —1.5 is the relevant index for ACDM 
around the non-linear scale, so we expect that for our universe the next-to-leading correction 
after the one-loop term is A 2 2 and not the two-loop correction as SPT would suggest. This 
should be checked directly by computing the EFToLSS parameters directly from simulations 
of our cosmology. 

6 A preliminary look at the n — — 2 plots, neglecting the unknown two-loop log term, shows that they are 
qualitatively similar to the n = —1.5, which we show in figure 2. 
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Figure 2: The plots show the fractional difference of EFToLSS (red, thick line) and, when available, 
SPT (blue, dashed line) predictions from the simulations of [8] for n = —1.5, —1 and of [9] for n = — 1. 
We show the same plot both with a linear and with a log abscissa. For n = — 1 when both simulations 
are available the top line refers to [9]. The /c-rangc shown and the ambiguity 0(1) in the definition of 
kjsiL are discussed in the middle of section 4. 

Let us move on and consider n = —1.5 and —1. In both cases there is one parameter, 
namely /3, which we fit as discussed above 7 . In the top of figure 2, we show n = — 1, for which 
the SPT prediction is cutoff dependent and hence unphysical. For this value of n both [8] 
and [9] provide fitting formulae so we show the fractional difference from both (the top red 
line refers to [9]). The EFToLSS provides a good fit in the dispalyed range. In the bottom 
panel, for n = —1.5 we compare with the simulations of [8]. In this case SPT is finite and 
so we show its prediction as a dashed blue line. One can see that the difference between 
SPT and EFToLSS increases as one approaches the non-linear scale as expected since the 
speed of sound and dissipation (and noise) terms grow with some positive power of k/k^L- 
The EFToLSS provides a better fit than SPT over the range shown. The results are not 

7 Notice that the parameter 7 comes from the noise term A^ oise . For the considered values of n, this is 
expected to be subleading to the two-loop correction, which we do not know. Hence it would be inconsistent 
to include it for n < —2/3. 
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yet conclusive because of the remaining uncertainty in the simulations, the small range of n 
studied etc. We postpone a more detailed numerical study to a future publication. 

5. Discussion 

Standard perturbation theory (SPT) is unsatisfactory for at least three reasons. First, it 
does not have a clear expansion parameter since the density contrast is not small on all 
scales. Second, it does not consistently account for deviations at large scales from a perfect 
pressureless fluid induced by short-scale non-linearities. Third, for generic initial conditions 
loop corrections are UV divergent and therefore the predictions beyond linear theory are 
cutoff dependent and hence unphysical. 

The recently developed Effective Field Theory of Large Scale Structure (EFToLSS) [2, 
3, 4] successfully addresses all of the above issues. The main idea is to smooth every field 
(e.g. density and velocity or momentum) by integrating out the short-scale degrees of freedom 
beyond a certain cutoff scale. If the cutoff scale is chosen to be larger than the non-linear scale, 
then the smoothed fields are small everywhere and provide a suitable expansion parameter 
for perturbation theory. The effect of integrating out the short-scale modes (i.e. substituting 
them with the solution of their equations of motion, in an explicit or implicit way) is to 
induce additional terms in the fluid equations for the large scales. At lowest order these 
are an effective pressure, shear and bulk viscosity and a stochastic noise, which shows that 
on large scales the behavior of collisionless dark matter deviates from that of a pressureless 
perfect fluid 8 . 

In this paper we have focused on the third issue mentioned above, i.e. the fact that SPT 
prediction are cutoff dependent and hence unphysical for a wide range of initial conditions. 
In the case of our universe, according to the ACDM model, these UV divergences are absent 
because the initial power spectrum is sufficiently red-tilted around the non-linear scale, with 
approximately n ~ —1.5. But this notwithstanding, one would like to describe our universe 
using a well-understood theory that is consistent with generic initial conditions and not just 
for specific ones. As our discussion here shows, investigating the theoretical consistency of 
the theory, even for phenomenologically non-viable initial conditions, highlights that we are 
missing some ingredient in the construction. In the present case, these ingredients are the 
effective coefficients induced by integrating out the short-scale. As pointed out before [2, 3, 4], 
and discussed here in more details, these terms (speed of sound, viscosity and stochastic noise) 
have precisely the right /c-dependence to cancel the divergences arising in perturbation theory 
beyond linear order. This implies that the EFToLSS, rather than SPT, is the consistent way 
to study the cosmological perturbation theory of structures. 

After we have convinced ourselves that these additional effective terms beyond a perfect 
fluid are needed for the consistency of the theory we can start asking about the role that 

8 These corrections, analogously to the non-linear corrections, go to zero as k — > 0. Neveetheless, if one 
wants to make accurate predictions, it is important to consistently include these corrections together with the 
non-linear ones, with the specific details depending on the initial power spectrum, as discussed in this work. 
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they play in phenomenological predictions. In order to assess the relative importance of 
various corrections we have studied here self-similar solutions in EdS. These are easy to handle 
analytically for a wide range of initial conditions and are not too far from realistic models. As 
we summarize in figure 1, depending on the index n in the self-similar initial power spectrum 
Pin = Ak n , the relative importance of the various corrections changes. For —3 < n < —2 
one- and two-loop corrections are the leading and next-to- leading order corrections 9 . For 
— 2 < n < — 1 one- loop is the leading correction and the speed of sound the next-to- leading 
one. This interval is expected to be the most relevant for phenomenology since in our universe 
around the non-linear scale n ~ —1.5. This is an important result since it shows that, for our 
universe, the corrections present in the EFToLSS, but neglected in SPT, are more important 
than two-loop corrections. A detailed study of the EFToLSS description of a realistic ACDM 
model was carried out in [3] (see also [4]). Here we have also compared the predictions 
of the EFToLSS with existing simulation of self-similar solution in EdS, in particular for 
n = —1.5 and —1. They fit the data well, but further analysis will be able to provide much 
more convincing evidence as we discuss below. As one goes towards larger n the EFToLSS 
corrections become even more important than the one- loop corrections. 

There are many interesting possibilities for future research. It would be appropriate to 
perform dedicated simulations of self-similar solutions in EdS to compare with the EFToLSS. 
Ideally, this would extend our results here towards larger n's for which the relative importance 
of EFToLSS increases. Also, one would like to focus on large scales where any subleading 
correction is very small, rather than close to the non-linear scale where all corrections tend to 
be of the same order. Other statistics beside the density power spectrum should be considered, 
and in particular comparisons performed realization by realization along the lines of [13]. One 
could study higher orders of perturbation theory with the goal of better defining the non-linear 
scale. 

Finally, in recent years it has become apparent that, at least for parameters close to those 
of our Universe, Lagrangian perturbation theory is a better starting point for a perturbative 
expansion [13]. The EFT approach has so far only been developed in the Eulerian context. 
It would be good to understand how the EFT ideas could be implemented in a Lagrangian 
scheme and thus perhaps improve it even further. 
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For n < —3 the theory is plagued by IR divergences. 
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A. Details of the dimensional regularization 



In [6] the results for the dimensional regularization of the P 22 and -P13 integrals were given 
for d = 3. From those formulae one can easily extract all the finite terms in the loop integral 
for those n for which there is no log-divergence. For the n's leading to a log-divergence in 
three-dimensions, we need to know how the loop integrals depend on d explicitely, in order to 
be able to take the limit e — > where e = d — 3 and extract the e _1 pole as well as the finite 
term e°. Hence in the following we re-compute the loop integrals in dimensional regularization 
and give and example of how the coefficient of the log-divergence is computed. 

A.l One-loop integrals 

In this section we re-derive the one-loop integrals following [6], but giving explicit formulae 
for the d dependence of the final result. For this computation we need the following formulae 
for the power spectrum (see e.g. [6]) with self-similar initial conditions Pi n (k) = Ak n 



P 22 (k,r) = A 2 a\r) 
P 13 {k,r) = A 2 a\r) 



d d q 

W) 

d d ci 

(2^ 



3 9 n 2|k 



qr[F 2 (s) (k-q,q) 



<f6ra (s) (k,q, 



(A.l) 
(A.2) 



with the kernels given by 
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where, instead of the general _F| (see e.g. [1]), we have written down the expression needed 
for the one-loop power spectrum. We will also need the mathematical identities 



/ 



d d q = rjd/2 - ^)T(d/2 - u 2 )Y{ Vl + v 2 - d/2) ^ d/2 kd _ 2ui _ 2u2 

(^2)^1 p _ ^2]^ 2 T(v x )T(v 2 )T(d -vi- v 2 ) 

k-q= ±J(fc 2 +g 2 -|k T q| 2 ) • (A.6) 



The computation can then be performed using Mathematica or a similar software. We attach 
to this publication a file [12] that outputs the final result. For P 22 one can straightforwardly 
use (A.6) with the upper sign to get P 22 in the form of a sum of terms like those appearing in 
(A. 5). For P13 one needs to use (A.6) with the lower sign for the terms with a factor |/c + gj -2 , 
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while (A. 6) with the upper sign for the others. The final result as function of d and n is 
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For d = 3 these general formulae reduce to 
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which is almost identical to the one in [6] except for an additional factor of (2tt)~ 3 , due to 
our different conventions for the Fourier transform, and the sign of the fourth term in P22, 
which we find to be minus instead of plus. Only with our sign the log-diverg ences in P22 arise 
for n = 1/2 + m, with m a positive integer, as expected from direct inspection of the cutoff 
integrals (see subsection 3.1). 

B. Self-similarity in EdS in d-dimensions 

In this appendix we collect some results about perturbation theory in d-dimensions in EdS 
(Q m = 1), which are useful when performing dimensional regularization. From the conser- 
vation of the matter energy-momentum tensor at zeroth order in perturbations one finds for 
the average energy density pa d = const and from the Friedman equation 

did ^H 2 = K( d y- P , (b.i) 

with K,(d) the d-dimensional Newton constant, one finds 

where for convenience we have set ciq = 1 and Tq -1 = 7io(d/2 — 1), so that one has % = 2/r 
in three dimensions. The non-linear system of fluid and Einstein equations before smoothing 
and neglecting vorticity, is given by [14] 

drSi + dt [(l + Stivf] =0, (B.3) 
d T v\ + Uv\ + d i( f> + vfd k vj = (B.4) 

didi(j)l = d i a 2 p5 l = d ^ 2 ±UH l , (B.5) 

i.e. only the Friedman equation is modified. As we did in section 2.2, one can immediately 
check that given a solution [5(x, t),v 1 (x, t), <f>(x, r)], another one is generated by the scaling 

S(x,t) = 5(X x x, \ t t) , (B.6) 
v i (x,T) = ^v i (X x x,X T T), (B.7) 

0(x,r) = (y) <t>(\ x x,\ T r), (B.8) 

for any A T and X x . The new solution corresponds to the same cosmology only if it has the 
same power spectrum of initial conditions. For this we need to know the growth factor. The 
linearized equation of motion for the density contrast is then given by 

-a 2 n 2 5, aa -^an\ a + d ^ J ^'H 2 5 = 0, or (B.9) 

-S trr - H5, r + ^^UH = , (B.10) 
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and therefore the two independent solutions are power laws: a d 2 ~ r 2 and a d / 2 ~ r d ^ 2 
In passing, the Green's function is then found to be 



G(a, a) = 6 (a — a) 



2d 



d-3 



H 2 (3d-4) 



ayd/2 / a \d-2 
a) \d 



(B.H) 



where 6(a — a) is the Heaviside theta function. The linear power spectrum, valid at early 
times and neglecting the decaying mode, depends on time as A 2 (/c,r) oc k d P(k,r) oc r 4 , 
i.e. the same in any dimension. The requirement that the rescaled solution represents the 
same cosmology as the original one is A 2 {k,r) = A 2 (k,r) which, given that A 2 (fc, r) = 
A 2 (k/X x , A r r), is satisfied if and only if 



\ x = \? 4+d . (B.12) 

for an initial power spectrum Pi n oc k n . So we conclude that in d-dimensions A 2 (k,r) has to 
be a function of k/k^i with 



NL ~ Aa 2 ( d - 2 ) A V r 
so that it is invariant under the rescaling (B.12). 
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